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ABSTRACT 

We develop a novel bosonic mean field theory to describe the spiral phases of a Heisenberg 
antiferromagnet on a one-dimensional chain, in terms of three bosons at each site. The ground 
state is disordered and for large values of the spin £, two different and exponentially small 
energy gaps are found. The spin-spin correlation function is computed and is shown to decay 
exponentially at large distances. Our mean field theory is also shown to be exact in a large- N 
generalization. 
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The study of the ground state and excitation spectrum of two dimensional quantum an- 
tiferromagnets has aroused considerable interest [1-14], particularly since their relevance to 
high T c superconductivity was realized [1,2]. A wide variety of approaches have been used 
- field theory methods [3-6], linear spin- wave theories [7], analogy with neutral superfluidity 
[8], bosonic and fermionic mean field theories [9-13], and numerical methods [14] - to study 
the problem. 

The first question that was tackled was the ground state and excitation spectrum of unfrus- 
trated Heisenberg antiferromagnets. Field theory methods showed that in one dimension, the 
ground state and its excitations had completely different properties depending upon whether 
the spin was an integer or half-integer. For integer spins, the ground state was exponen- 
tially disordered at long distances and excitations had a gap, whereas for half-integer spins, 
the ground state was only algebraically disordered and had massless excitations [3,4]. This 
difference had its origin in the topological term that was induced in the long wave-length 
effective field theory. But this difference disappeared in two dimensions, where no topological 
term was found to be induced [5]. Moreover, numerical evidence [15] favoured a Neel ordered 
ground state in two dimensions. 

The question of the ground state and excitation spectrum of a frustrated model is far more 
complex and is, as yet, not completely understood. The mapping to conformal field theories 
[4] has led to the expectation that the phase diagram for half-integer spins has a region 
where the model has algebraic disorder and massless excitations. Outside this region, the 
system is expected to be dimerized, except at specific points. For integer spins, the region 
of massless excitations is replaced by massive excitations. However, there exist few explicit 
results. For two dimensional models, the phase diagram is even more uncertain, because the 
various different approaches lead to different answers. In particular, the existence of a spin- 
liquid state (predicted by fermionic mean field theories) is still unconfirmed. (For a recent 
summary of existing results, see Ref. [16].) 

In this paper, we introduce a new (bosonic) mean field theory (MFT) involving a repre- 
sentation of spins by three bosons in the adjoint representation of the SO (3) group of spins, 
precisely to study this question. In particular, we address the following specific issue. We 
study the frustrated Heisenberg antiferromagnet (AFM) on a chain using our three boson 
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representation (3BR), with the aim of obtaining explicit results. This method (like all mean 
field methods) is insensitive to the presence of topological terms. However, since the ultimate 
aim is to study spirals in two dimensions, where topological terms are not expected, this is 
not a serious handicap. This 3BR works well for spiral phases, reproducing the zero modes 
at q = and q = 9 (where 9 is the spiral angle) as expected from symmetry considerations. 
The same representation also works for the spiral phases of the triangular AFM in two di- 
mensions and is expected to work for the helicoidal phases of the frustrated model on a two 
dimensional square lattice. One advantage of the 3BR is that no rotation of all the spins to 
a ferromagnetic configuration is needed. This makes an analysis of helicoidal phases in two 
dimensions technically much simpler [17]. 

In the 3BR employed here, we represent the spins at every site by a triplet of bosons. To 
enforce the spin nature of the operators, two constraints are required at each site which is in 
contrast to the more commonly used representation of spins in terms of two Schwinger bosons 
(2BR) where only one constraint per site is required. These two constraints are imposed on 
an average. Following the method of Sarker et al[12], we perform a Hartree-Fock-Mean-Field 
(HF-MF) averaging to obtain the ground state energy in terms of six variational parameters, 
whose values, in turn, are obtained by extremizing the energy. In the S — > oo limit, the spin- 
wave spectrum is reproduced. For S large but finite, our solution yields two exponentially 
small energy gaps. We can compute the two spin correlation function and show that beyond 
some length scale, the correlation function falls off exponentially. Finally, we show that our 
HF-MF treatment becomes exact in the large- N limit by generalizing the 3BR (forming the 
triplet representation of SO (3)) to an iVBR (forming the iV-plet representation of SO(N)). 

Let us start with the frustrated Heisenberg AFM on a chain described by the Hamiltonian 



where we have normalized the exchange constant J = 1. For 6 < 1/4, the classical ground 
state is Neel ordered, whereas for 5 > 1/4, the classical ground state is a spiral where all the 
spins lie in a plane and the relative angle 9 between any two spins is given by 




(1) 



n 



cos 9 = -1/45 



(2) 
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so that 7r/2 < 9 < tv. The classical ground state energy per unit spin is given by 

to leading order in S, where N is the total number of spins. For 5 = 1/2, (a special case 
called the Majumdar- Ghosh model [18] ), this model was recently studied by us [19]. In the 
long distance, large-*? limit, we mapped the model to an .SO (3) -valued field theory and using 
the /3-functions of the theory, we showed that the ground state was exponentially disordered 
and exhibited a gap. However, we were unable to generalize that method to arbitrary 5. (See 
also Ref. [20] ). 

We first perform a spin- wave analysis of this Hamiltonian using Villain's action-angle vari- 
ables [21]. For the general spiral case - i.e., 5 > 1/4, - the spin-wave spectrum [22] (valid for 
large- S) is given by 

ujk = 2S [ ( — cos# — 5 cos26> + coska + 5 cos2/ca) 

(4) 

( — cos6> — 6 cos26> + cos6> coska + 6 cos26> cos2/ca) ] 1 ^ 2 . 

(The lattice spacing is a). Within the first Brillouin zone, u>k vanishes at ka = 0, +9 and —9 
with a linear dispersion, and the spin-wave velocities are given by 

c = Sa (1+45) (1-1/165 2 ) 1/2 

(5) 

and c e = c (1 - 1/25 + 1/85 2 ) 1/2 

at ka = and ka = ±9 respectively. Thus, the mode at k = has a higher velocity than the 
two modes at ka = ±9. We shall see that this spectrum is reproduced by our bosonic MFT 
in the large- S limit. 

We now set up the bosonic MFT by expressing the spin in terms of bosons. The 3BR 
expresses the components of a spin S in terms of three bosons as 

S a = -i e a/ 3 7 aj; a 7 (6) 

where a, (3 and 7 run from 1 to 3, e Q/ 3 7 is completely antisymmetric and repeated indices are 
always summed over. Using [a Q ,a^ ] = 5 a p, we can check that the spin algebra is satisfied. 
But to also have S a S a = S 2 = S(S + 1), we need to impose the constraints 

a a a a = S and a a a a ap ap = (7) 
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on any physical state. Notice that the first equation in Eq. (7) implies that the 3BR works 
only for integer spins. By enforcing these constraints, we can check that for any S, the total 
number of orthonormal states is 2S + 1 as expected. The second constraint in Eq. (6) can 
equivalently be rephrased as the two constraints 

a a a a = an d = (8) 

in the sense of matrix elements between any two physical states. This is the form in which it 
is employed later. 

To understand the connection between spin order parameters and appropriate expectation 
values of the bosons, consider a spin operator lying in the x-y plane, in the S — > oo limit, - 
i.e., in the classical limit, - with expectation value 

< S a > = S (cos0, sin0, 0). (9) 

Then Eqs. (5) and (6) imply that 

< ai > = i sin0, < a 2 > = -i cos0, < a 3 > = (10) 

upto an arbitrary phase multiplying all the bosons. So if (f) mn is the angle between the m th 
and n th spins, we have 

< S m • S n > ~ < S m > • < S n > = S 2 cosc/) mn (11) 

which in turn implies that 

< a ] ma a na > ~ < a ] moi > < a nce > = S cos 2 (0 mn /2) 

(12) 

< a ma a na > ~ < a ma > < a nce > = S sin 2 (0 mn /2), 

again upto arbitrary overall phases. Thus, spiral ordering of spins implies non-zero expec- 
tation values of bosonic bilinears. In fact, < a^ana > is the Ferromagnetic (FM) order 
parameter and < a ma a na > is the Antiferromagnetic (AFM) order parameter. 
Next, we observe that the product S m • S n can be written as 

S m " S n = : Ym,n Ym,n '■ — -^-m,n -^-m,n (13) 



where 

Y m ^ n ^rna ^nct and X mn (t ma ^na- (^^) 

In terms of the bilinear s X m ^ n and Y m ^ n , the Hamiltonian in Eq. (1) can be rewritten as 

n 

+ A n (a^Q, a nQ , — S) — p n (a. na a na ) — p n a^ a ], 

(15) 

with A n ,p n and p* being the Lagrange multiplier fields introduced to enforce the constraints 
in Eqs. (7) and (8) at each site. We now make a HF decomposition by writing 

A* A = < > A + A* <A> - < A ] >< A >, (16) 

where A = X n>n +i,Y n>n+ i,X njn+ 2 and Y n ,n+i in turn. Such a decomposition (in contrast 
to the Peierls variational decomposition, which allows for all possible decouplings where each 
four boson term is written as products of pairs in three different ways) is justified later 
by a large-iV argument. Further, we make the MF ansatz that the parameters A n = A, 

Pn = Pn = < X n,n+1 > = X\ , < X n ^ n+2 > = X 2 , < ^n,n+l > = Y\ and < F n>n+2 > = 

Y 2 are all independent of n and real. Notice that this ansatz - in particular, the reality of p, 
Xi and Yi - breaks the local gauge invar iance a na — > e %dn a na of the Hamiltonian. However, 
physical quantities such as spin-spin correlations remain gauge-invariant [12]. Also note that 
non-zero values for Xi, X 2 and Yi, Y 2 imply the existence of short-range AFM and FM 
orderings respectively. 

We now diagonalize the MF Hamiltonian by a Bogoliubov transformation to obtain 



Hmf -= _ XS + X\ + 8Xl - Yi - SY 2 2 



N 



r /a dk , t t l 3 3 n 

+ a — ( u k b k b ka + - u k - - fj, k ) 

J n 2 2 



(17) 



where 



/ 2 2 \l/2 

= ( Pk ~ "k ) ' 

H k = A + 2Yi cosfca + 25 Y 2 cos2/ca (18) 
and v k = p + 2Xi coska + 25 X 2 cos2/ca, 
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and the bosons b kct are related to the bosons a ka by the standard Bogoliubov transformation 
given by 

a ka = cosh6> fc b ka + sinh6> fc b^_ kce 



a\ a = sinh0 fe b- ka + cosh6> fc b\ a . 



(19) 



The factors of three in Eq. (17) arise because the constraints in Eq. (7) have only been 
imposed on the average, resulting in the decoupling of the three bosons. (This tripling of 
number of branches of the spin-wave spectrum is an unfortunate feature of this MFT. A 
similar doubling of number of branches occured in the MFT based on the 2BR [11] [12].) 
Thus, we obtain the MF ground state energy as 

= - (S + 3/2) A + X\ + 5X1 - Y? - 5Y? + ^ f * u k . (20) 

The equations for the six variational parameters A, p, Xi and Yi are obtained by extremizing 
the energy and are given by 

s+- - — r /a — ?± 

2 2 Jo tv u k 
Q = 3a r' a dk VjL 

2 Jo 7T UJ k 

Y - ^ ^ d - !* cos ka 

2 Jo * (21) 

v 3a r' a dk v k 

X\ = — / — — cos ka 

2 Jo 71 u k 

3a r' a dk ii k 
Y 2 = — / — — cos 2ka 

2 Jo 71 u k 

, Y 3a r' a dk v k 
and A 2 = — / — — cos 2/ca. 

2 Jo 71 u k 

These equations look rather intractable. But, in fact, it is possible to obtain the solutions 
to leading order in S. We know that as S — > oo, the solution should approach the classical 
spiral ground state configuration asymptotically. Hence, to leading order in S, we must have 

Yi = S cos 2 0/2, X 1 = S sin 2 6/2 

Y 2 = S cos 2 9, X 1 = S sin 2 (22) 
A = - 2S (cos# + 5 cos 29 ) and p = 0, 
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where 9 is defined in Eq. (2). (One can check that these are the correct asymptotic values of 
the parameters by comparing them with the values given in Eq. (12) as well as by comparing 
the dispersions given in Eqs. (18) and (4). Note also that by substituting this solution in Eq. 
(20), the MF energy agrees with the classical energy to leading order in S.) However, with 
these asymptotic solutions, notice that u)k — > at k = and k = 9/a. Hence, the right hand 
sides (R.H.S.) of Eq. (21) are log divergent and can only equal the left hand sides (L.H.S.), 
which are large but finite, if we allow for a small mass generation. Let us assume that near 
k ~ and k ~ 9, the dispersions of ui^ are given by 



oo k ~ \ A 2 + elk 2 , k ~ 

r-L (23) 

and Uk — yAg + c 8 (k — 9/a) 2 , k ~ 9/a, 

where cq and eg are the spin- wave velocities given in Eq. (5), and Ao and Aq are the two 
small masses generated. In fact, we will show that the A; are exponentially small - i.e., of 
0(e~ s ). Notice also that the non-singular regions on the R.H.S. of Eq. (21) are of 0(1) 
(not of O(S) ) and do not contribute to establishing the equality of the L.H.S. and R.H.S to 
O(S). Hence, to leading order in S, we can simply assume that the integrals on the R.H.S. 
are dominated by their values at k = and k = 9/a. In this limit, Eqs. (21) reduce to 

S + 3/2 = (J>o I + He/a Ig/ a 
= Vq Iq + VQ/a h/a 

Yi = no I + He I a cos/9 I e/a 

Xt = i/ I + v e /a cos6> I e/a 

Y 2 = Ho J + He/a cos20 I e/a 

and X 2 = u Q I Q + u 9/a cos26> I e/a , 



(24) 



with 



3a f e dk 1 3a eco 

Jo = — / — — — m — — 

2 J n u k 2-kco A 

, T 3a f / a+£ / 2 dk 1 3a , ec 

and 10/ a = — — m ——. 

2 J 0/a _ e / 2 7T 27TC e A e 



(25) 



where A/c << e << n/a . We can now explicitly check that the O(S) terms on both sides 
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of Eq. (24) are satisfied when 



c r 2tvS 45-1 

/\ Q g X p 



and Ag 




a ^ L 3 \ 45 + 1 1 

7T-, (26) 

45-1 



(45 + l)(l-l/25 + l/85 2 ) 



To summarize, the solutions to the Eqs. (21) are given in Eqs. (22) and they lead to two 
exponentially small mass gaps in the theory given in Eqs. (26). 

It is interesting to compare the values for the gaps in Eqs. (26) for 5 — 1/2 with the values 
obtained for the same 5 using the field theory approach [19]. The one-loop /3-function of 
the field theory of the Majumdar-Ghosh model led to the single mass gap A ~ exp(— l^S 1 ), 
whereas here, the bosonic MF treatment yields A ~ exp(— 1.2S) and ~ exp(— .865 1 ), 
which agree, at least upto the order of magnitude of coefficient of S. The two mass gaps 
obtained here (like the two spin-waves) appear to reflect the fact that the spiral ordering of 
the ground state picks a particular plane. Thus, it is not surprising that fluctuations, and 
hence, onset of disorder, within the plane and perpendicular to the plane, have different mass 
scales. The field theory method, presumably, was not sensitive enough to see this feature. 

Let us now calculate the spin-spin correlation function within the bosonic MFT. From 
Eqs. (13) and (14), we see that the product of any two spins can be written as a product of 
bilinears, so that the spin-spin correlation function is given by 

< So • S n > = <: Yl n Y , n :> - < X{ n X , n > . (27) 

Using only the Wick contractions allowed by the HF decomposition in Eq. (16), we find that 
the spin-spin correlation can be written as 

< So ■ S n > = | < Yb,n > | 2 — I < -Xo,n > | 2 , (28) 

where < Yb,n > and < X$ >n > are obtained by using the Bogoliubov transformation in Eq. 
(19) as 



3a r /a dk ,Hk 
< Y , n > = — / — ( — - 1) cos nka 
* Jo n uu 



(29) 

3a r /a dk v k 



. r ■' dk „ 

and < X n > = — / — — cos nka 

2 Jo K uj k 
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Once again, the integrals are dominated by the regions near k ~ and k ~ 9. We now 
explicitly compute the correlation function in two limiting cases. When na, the distance 
between the two spins measured in terms of the lattice spacing a, is small, - i.e., na << 
Aq 1 , Ag 1 , we find that Y"o,n ~ Scos 2 (nO/2) and X , n ~ sin 2 (n0/2) so that 

< S • S n > = S 2 cos n6, for na « Aq 1 , A^ 1 . (30) 

This is not surprising, because at short distances, we expect the system to be ordered. But 
at long distances, - i.e., when na » Aq 1 , A^" 1 , - we find that 



n [ e dk cos nka _ T)0 A n /rn 
<y ,n> + <X , n > ~ 5 — - ~ Se naA °/ C0 ,and 

Jo ^ v /"fc — ^fc 

^ rr f e/a+e/2 dk cos nka _ nflAll / c „ 

< y ,n > - < X ,n > ~ S / — . : ~ 5 cos n# e naA "/ ce 

Je/a-e/2 K V^fc + ^ 



(31) 



Hence, for large enough distances, 

< So • S n > ~ S 2 cosnO exp [—na( — - H -)], for na >> Aq 1 , A^" 1 , (32) 

Co eg 

- i.e, the correlation function falls off exponentially. 

Let us now justify the HF decompositions or Wick contractions used in Eqs. (16) and (27). 
Naively, the four boson terms that appear in the product of two spins can be decomposed (or 
contracted) in three different ways. However, in our MF treatment, we have only allowed one 
possible decomposition. For example, n X miTl = al na a^ a a m/ 3 a n p is only decomposed as 
< a lna. a na > a mpCL n /3 + ^lna a lia < a m(3a n /3 >■ The other possible contractions < a^ a nf3 > 
and < a\ na a m p > which are down by factors of 1/3, because the 50(3) indices a and (3 are not 
summed over, have been ignored. The justification for this treatment comes from the large- iV 
generalization of the model. The 3BR is generalized to an iVBR (N boson representation) of 
SO(N) 'spins'. In the N — > oo limit, the other contractions which are down by a factor of 
1/N can certainly be ignored. 

Our large-iV generalization is similar in spirit to the generalization of the 2BR of 577(2) 
spins to the N boson representation of SU(N) spins discussed in Ref. [11]. But just for 
completeness, we mention some details of our large- TV model. We write the components of 
an SO(N) spin S a as 

S<* = -i 9^ <4 a 7 ( 33 ) 
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with a = 1, . . . N(N — l)/2 and (3, 7 = 1, . . . N. Furthermore, we have the relations 

9pj = - 9jp and gf € = 5p 5 8 le - 5p t 8 l5 . (34) 

To reproduce the spin algebra and the correct number of states, the constraints in Eq. (18) 
are now replaced by 

UaUot = NS/3 and a^a^apap = 0. (35) 

Clearly, this representation works only if NS/ 3 is an integer. The Hamiltonian for general N 
given by 

3 . 

H = — ( S n • S n+ i + 5 S n ■ S n+2 ) (36) 

n 

can be written in terms of the bosons in Eq. (33) along with the Lagrange multiplier fields to 
enforce the constraints just as was done earlier. In fact, the entire analysis can be reproduced. 
Here, however, our aim in introducing the large-iV formalism was only to justify the HF 
decoupling procedure that we used. 

To conclude, a notable feature of our analysis is that we do not need to rotate the spins of 
the ground state of interest in order to make it look FM and then proceed with the analysis. 
For the 2BR, such a rotation is usually performed [11]-[13]. We are currently using the 3BR 
to study frustrated spin models in two dimensions [17]. 
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